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Abstract 

We consider a planar Brownian motion starting from O at time t = and stopped at 
t = 1 and a set F — {Oli ; i = 1, 2, . . . , n} of n semi-infinite straight lines emanating from 
O. Denoting by g the last time when F is reached by the Brownian motion, we compute the 
probability law of g. In particular, we show that, for a symmetric F and even n values, this 
law can be expressed as a sum of arcsin or (arcsin) 2 functions. The original result of Levy 
is recovered as the special case n — 2. A relation with the problem of reaction-diffusion of a 
set of three particles in one dimension is discussed. 



The first arc-sine law gives the distribution of the number of positive partial sums in a 
sequence of independent and identically distributed random variables. It was first discovered by 
P. Levy in his study of the linear Brownian motion and then discussed a lot for its relevance to 
the coin-tossing game 1 . The second arc-sine law, also discovered by P. Levy [2], provides an 
information on the last passage time which can be stated as follows. Consider a linear Brownian 
motion B(t) starting at at time r = and stopped at time t and let g be the last time when 
is visited. The random variable 

5 = sup{r <t,B(r) = 0} (1) 

satisfies 

2 fu . . 

P(g < u) = — arcsin J— (2) 

with the density V(u) = = (3) 

vr yu(t — u) 



Over the years, this result has been extended in several different directions (see for instance, 
[31 E]) and is still a subject of active research in probability Generalizations of the first arc- 
sine law have also been considered in different contexts (one-dimensional diffusion in a random 
medium 6 , Brownian motion on graphs [Zj and, also, in two dimensions |S]). 
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The purpose of this Letter is to present a two dimensional generalization of the law (|2"j). As a by- 
product of this result we also derive an explicit expression of the first passage time distribution 
which is relevant for a problem of reaction-diffusion involving three identical particles. 

Exit problems for Brownian motion have a rich history and several applications in physics 
(see for instance They are in particular related with problems of capture of independent 
Brownian particles diffusing on the line. This connection was first anticipated by Arratia jlUj 
and then discussed in the mathematics I12j and physics litterature I13[ I14| I15j mainly in 
the context of reaction-diffusion models. In the case of three particles, the process (x\(t) — X2(t), 
X2(t) — xs(t)) defines a certain diffusion in a quadrant of R 2 . By a suitable transformation, this 
process can be mapped on a diffusion inside a wedge whose angle depends on the diffusion 
constants. Using this correspondance, it has been shown that the first passage time through 
the wedge gives the survival probability; a quantity which decays with a power law which only 
depends on the angle of the wedge |14l I16| . At the end of this work, we exploit this correspon- 
dance to compute exactly (and not only asymptotically) the first collision time distribution for 
a three particle problem. Our approach is based on an identity relating the first passage and 
last passage distribution which has an interesting probabilistic interpretation [T7] . 

To begin with, let us start by considering, as in Figure 1, a wedge of apex O and angle (ft 
with a boundary F = {OJ, ; i = 1,2} and a two dimensional Brownian motion ~r*(t) starting 
from O at t = and stopped at t somewhere in the plane. We denote by g the last time when F 
is visited and compute the probability P(g < u). Due to the scaling property of the Brownian 
motion, this distribution is a function of the reduced variable u/t. In the following, we will for 
simplicity set t = 1. 
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Figure 1: The frontier F = {0/1,0/2} of the wedge divides the plane into 2 regions, (1) and 
(2). 



Suppose that the particle reaches some point Tq at time t = u (see Figure 1). Clearly, if Tq 
belongs to region (1) (resp.(2)), the particle must stay in (1) (resp.(2)) between t = u and 
t = 1 in order to satisfy the condition g < u. We can therefore write 

P{g<u) =pW(u)+pM(u) (4) 

Expressing the fact that the propagation is free between t = and t = u and that the particle 
has not hit the boundary between t = u and t = 1, we get 

P®(u)= [ d 2 n ( d 2 T — G®(T,1; r$,u) ; i = l,2 (5) 

J(i) J(i) 2-7TU 
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The propagator satisfying the diffusion equation with Dirichlet boundary conditions on F 
is given by 



Ml) 2 ^ ■ m07T ■ me o 7r -£tH t { rr o \ , R \ 
G ( ' = — — r > sin — — sm — ; — e I mm. (6) 

0(i-«) * u-«/ 

where I u is a modified Bessel function and the notations are defined on Figure 1. 
Performing the spatial integrations in (JSJ), we get [TH 

P«r ^ 1 v [ r (^ + fc + rf i m 



P^ 2 ) is obtained by the change eft — > 2tt — (ft in (JJJ). Therefore, for arbitrary values of the 
law P(g < u) is written in terms of a double series. 

As a check, let us first consider the special case eft = n. We may write 

p,fc=0 v ' 

2 ._ 

= — arcsinvu, (9) 

7T 

The fact that one recovers Levy's second arc-sine law is not surprising since, when (j) = n, F 
divides the plane into two half-planes. Therefore the component of the Brownian motion parallel 
to F factorises and plays no role: we are thus left with a one dimensional problem. 



Coming back to general values of </>, we can derive the behavior of the probability density V 

du 



dP(g<u) ^ w k en u _^ q anc j u _^ ^ gy us i n g Q anc [ (|7|). one gets a power-law behavior when 



0+ 



T(u) ~ - n~ 1/2 for = vr (10) 

7>(u) ~ C(/i) nt" 1 for 4> ^ 7T (11) 

with c <"> = ^i^TI) (12) 

and /U = when < <p < n (13) 

2-7T — cp 

fi = — when 7r < < 27r (14) 



Now, for the limit u — > 1 , using asymptotic expansions for T functions and also an equivalence 
between series and integrals, we get 

Viu) ~ - — L= (15) 

i.e. the same behavior as for ©. The expression (|15|) doesn't depend on <p and we have already 
seen that 4> = tt gives the Levy's law. Remark that u —* 1~ corresponds to Brownian curves 



3 



that stop close to F. Therefore, between t = u and t = 1, the Brownian particle only "sees" an 
infinitesimal part of F, i.e. a straight line as for (ft = ir. This is, in our opinion, why the result 
(|15|) doesn't depend on (ft. Actually, it only depends on the fact that the plane is divided by F 
into two regions. We will come back to this point latter on. 




Figure 2: F consists in n semi- infinite straight lines starting from O. 



To go further, let us remark that for F = {Oli ; i = 1,2, . . . , n} as in Figure 2, Q) becomes 
simply: 

n 

P(g<u) = Y J P {€) {u) (16) 

i=l 

(Replace (ft by (fti in (JJJ) in order to get P"). 

Let us now specialize to the situation when F is symmetric and n is even (n = 21). In that case, 
F consists in I infinite straight lines crossing at point O and dividing the plane into 21 equal 
angular sectors, each one of angle (ft = it /I. 

a) 1=2 b) 1=3 



O 




Figure 3: F is symmetric. The analytic form of P(g < u) will depend on the parity of I. Thus, 
it will be different for the cases a) and b). For further explanations, see text. 
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Equation (jlflj) writes 



p f„ fcS (2ij- + ! + *)! *l 



Z being an integer, we can sum the series and, finally, get 



fc=o 



fc=0 



(17) 



21 ( 2irk\\ 

Pl(u) = — arcsin ^-v/ucos — ^— J J , /odd (18) 

P^u) = ^ ( _1 ) fc (^arcsin ^^/ucos^jj j , Z even (19) 
which is the central result of this paper. 

We remark that the correct small u behavior for Pi(u) follows from the two identities 
l ~ l / 27rk\ m 

^(cos — ) =0 , Z odd , m = 1,3, 2 (20) 



J2(-l) k (cos^) = , I even , m = 0, 2, 4, 2 (21) 



In particular 



P A (u 
P 5 (u 



arcsm v u 



IT 

4 



arcsm 



V^) 2 = Pf 



u 



6 / ^- . v , , 

— arcsm y u — 2 arcsm 

7T \ 2 



7T" 



(arcsiny^) 2 — 2 ^arcsin y ^ 



10 / _ / 7T _\ / 2vr _ 

— I arcsin — 2 arcsin I cos — yu\ + 2 arcsin I cos — yju 



These functions are displayed below in Figure 4. 



(22) 
(23) 

(24) 

(25) 
(26) 



As expected, the Levy's second arc-sine law is recovered in ()22|) , Moreover, the result (|23|) is 
straightforward since, when I = 2, the 2 components of the Brownian motion factorize. Thus, 
for I = 2, (arcsin) 2 functions appear. What is surprising is that they will appear each time I is 
even while being absent when I is odd. 



For the probability density, V\ (=757-), with (fTH |) and ([T9Jl . we obtain: 



Viiu) 



I 



vr y/l - u 

This is consistent with (|15jl that corresponds to / = 1 



when 



u -» 1" 



(27) 
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We now present a formula which relates the first passage and the last passage time distribu- 
tion. The starting point is Q) and © which may be rewritten as 

f 1 r ^ 

P( g <u)= Pr(T > 1 - u\r ) - r dr (28) 
J u 

where Pr(T > (1 — u)\ro) is the probability distribution of the first passage time T through F, 
given that the process starts at rQ. Then, by scaling one has 

Pr(T > (1 - u)\r Q ) = Pr{T > ^ ~ ^ |1) (29) 

By a simple change of variables it follows that 

p ^ < = J Pr(T > e ~ x dx m 

Therefore 

P(9 < Y^~ t ) = E(e-J?) (31) 

which is a relation between the first passage characteristic function for a process starting at 
ro = 1 and the probability distribution of the last passage time. Interestingly enough this 
formula can also be derived in a more intrinsic fashion using only time inversion and scaling 
[T7| . As an application, let us the derive the density of first passage time in a wedge of angle 
(f> = ^. In the context of the capture problem mentioned in the introduction, this corresponds 
to a set of three identical and independent particles In this case, the distribution P(g) is 
given in eq (J24[). By an inverse Laplace transform (|31l) gives the density of first passage time: 

f(T) = 4-e-*( ey'dy-2 e^dy) (32) 

TTlT JO JO 

One can check that this formula is in agreement with (16) of which expresses the first 
collision time probability for a given set of initial conditions. By averaging this formula over the 
angle and setting r = 1 one recovers ([32 [I. 
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